Logarithmic Relaxation in a Kinetically Constrained Model 
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Very recently, we have reported a numerical inves- 
tigation on a simple bead-spring model for polymer 
blends with components of very different relaxation 
timesi. At moderate and high dilution of the fast 
component, the slow component acts as a confining 
medium for the former. In this situation, dynamic cor- 
relators for the fast component display non standard 
relaxation features. Instead of the usual two-step de- 
cay observed for dense colloidal systems or deeply su- 
percooled liquids, dynamic correlators do not exhibit 
a defined plateau and the long-time decay is highly 
stretched. Decreasing temperature yields a striking 
concave-to-convex crossover in the shape of the corre- 
lators. At some intermediate temperature, the decay 
is purely logarithmic over 3-4 time decades. Atom- 
istic simulations on real polymer blends show a quali- 
tatively similar behavior—. These results exhibit strik- 
ing analogies with predictions of the Mode Coupling 
Theory (MCT) for state points close to a higher-order 
transition'^''*. In particular, asymptotic laws predicted 
by this MCT scenario for dynamic correlators provide 
an excellent description of simulation results for the 
fast component. 

Higher-order MCT transitions have been strictly 
derived for simple models of short-ranged attractive 
colloids^i^iSii. For these systems, such transitions arise 
in regions of the control parameter space where re- 
versible bond formation (originating from the short- 
ranged attraction) and hard-sphere repulsion compete 
as different mechanisms for dynamic arrest. The ob- 
served strong dynamic analogies with polymer blendsi 
suggest that this anomalous relaxation scenario might 
be a general feature of systems where dynamic arrest 
is driven by several competing mechanisms of different 
characteristic lenghts. For the fast component in poly- 
mer blends we have suggested confinement and bulk- 
like cagingi-, respectively induced by the host matrix 
formed by the slow component, and by neighboring 
particles of the fast component. 

Kinetically constrained models^-^ are used as 
coarse-grained pictures for relaxation in supercooled 
liquids. In these models, the liquid structure is sub- 
stituted by an array of cells. The cell size roughly 
corresponds to a density correlation length. A spin 
variable is assigned to each cell, with values or 
1 denoting respectively unexcited and excited local 
states in a mobility field. Change in local mobility 
(spin flip) for a given cell is permitted according to 
kinetic constraints determined by the mobilities of 
neighboring cells. Propagation of mobility, leading 
to structural relaxation, occurs via dynamic facilita- 
tion: microscopic regions become temporarily mobile 
only if neighboring regions are mobile. Kinetically 
constrained models provide in a simple way an impor- 
tant feature of glass-forming liquids: the increasing 
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FIG. 1: Temperature dependence of the mean persistence 
function P{t) at composition X{ = 0.33. Symbols cor- 
respond to fast cells. Dashed lines correspond to slow 
cells (same temperatures as for the former ones, decreas- 
ing from left to right). Solid straight lines indicate nearly 
logarithmic relaxation. 



of dynamical correlation lengths with decreasing tem- 
perature, leading to dynamic heterogeneitjiifl*ii*iSiii. 

Some of the anomalous dynamic features displayed 
by short-ranged attractive colloids, as logarithmic re- 
laxation or reentrant behavior of the diffusivity, have 
been recently reproduced by a simple model based on 
this latter approach^. Motivated by this fact and 
by the mentioned dynamic analogies between short- 
ranged attractive colloids and polymer blends, in this 
Note we investigate a simple kinetically constrained 
model aimed to reproduce dynamic features in real 
binary systems with components of very different re- 
laxation rates. Results reported here exhibit strong 
analogies with dynamic correlators in polymer blends. 

We have carried out Monte Carlo (MC) simulations 
on a variation (see below) of the north-or-east-or- front 
(NEF) model recently investigated by Berthier and 
Garrahan^^. In these three-dimensional model, spin 
flip of a given cell is only permitted if there is at least 
one excited neighboring cell in the north-or-east-or- 
front directioni^. More specifically, if the cell is de- 
noted by its position vector i = (ix,Jy,Jz), spin flip is 
permitted according to the following rules: 

i) At least one of the neighboring cells in the north 
(ix,iy,Jz + l), or east {ix,iy + l,iz), or front {i^ + l,iy,iz) 
direction is excited, i.e., it has spin 1. 

ii) If i) is fulfilled, spin flip is accepted according 
to the Metropolis rulei^. Hence, 1 ^ is always ac- 
cepted, while — > 1 is accepted with a thermal prob- 
ability [1 + exp(l/T)]~^, where T is the temperature. 

In the present work we investigate a binary mix- 
ture of cells with the same population of excita- 
tions, [1 -|-exp(l/T)l~*, but different rates for excita- 
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FIG. 2: Symbols: composition dependence of P{t), at T = 
0.50, for tiie fast cells. The straight line indicates nearly 
logarithmic relaxation. 



tion, exp{—E/T)[l + exp(l/T)]^"'^, and de-excitation, 
exp(— S/T), of cell mobility. This choice fulfills de- 
tailed balance. We use E = and 3 for the diff'erent 
components, which are respectively denoted as "fast" 
and "slow" cells. Hence, spin flip rules for the fast 
component are the same as in the original NEF model, 
aimed to reproduce bulk-like dynamics. The probabil- 
ity of spin flip for slow cells is decreased by a factor 
exp(— 3/r) as compared to that of fast cells, providing 
a large separation in the time scales of both compo- 
nents (see below). 

The mixture composition, X{, is defined as the frac- 
tion of fast cells. We investigate a wide range of com- 
position and temperature as control parameters. A 
square box of side A^ = 50 cells is used. Periodic 
boundary conditions are implemented. Slow and fast 
cells are randomly distributed according to the se- 
lected value of Xf . Time is given in MC steps. Within 
each MC step a total of N^ trials (one for each cell) 
are performed. 



We have evaluated the mean persistence function 
P{t) = "^iPii; t), where p{i; t) is the persistence func- 
tion of the cell i at time t. It takes the value 1 if 
no spin flip has occurred for that cell in the interval 
[0,i]. If one or more spin flips have occurred, it takes 
the value 0. Fig. ^ shows results for P{t) at com- 
position Xf = 0.33 and different temperatures. Note 
that the first decay usually observed for any correla- 
tor in real systems — corresponding to the transient 
regime at microscopic times — is obviously absent due 
to the coarse-grained character of the model. The 
introduction of very different rates for cell spin fiip 
provides very different relaxation times for fast and 
slow cells. Slow cells display standard relaxation, as 
observed for the slow component in polymer blends^'. 
However, fast cells display rather different relaxation 
features. In complete analogy with results for the 
fast component in polymer blendsi^, P{t) exhibits a 
concave-to-convex crossover with decreasing tempera- 
ture. At an intermediate temperature, pure logarith- 
mic relaxation is observed over a time interval of near 
three decades. This behavior is also observed by fix- 
ing temperature and varying the mixture composition 
(see Fig. (SJ. It is worthy of remark that the features 
displayed in Figs. ^ and [21 are not specific of P{t). 
Qualitatively analogous results are exhibited by other 
dynamic correlators, as persistence correlators for dis- 
tinct cells, or intermediate scattering functionsii. 

In summary, the highly non-trivial anomalous re- 
laxation features recently reported for the fast com- 
ponent in polymer blends can be qualitatively repro- 
duced by a simple kinetically constrained model. The 
fact that another recent work within this approachi^ 
also reproduces unusual relaxation features for short- 
ranged attractive colloids suggest the relevance of the 
dynamic facilitation picture for systems with several 
competing mechanisms for dynamic arrest. 

We thank J. P. Garrahan and L. Berthier for useful 
discussions. 
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